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Abstract 

Background: To examine the effects on LDL cholesterol of the combined use of statins and phytosterolsAstanols, 
in vivo studies and clinical trials are necessary. However, for a better interpretation of the experimental data as well 
as to possibly predict cholesterol levels given a certain dosing regimen of statins and phytosterolsZ-stanols a more 
theoretically based approach is helpful. This study aims to construct a mathematical model to simulate reductions 
in low-density lipoprotein (LDL) cholesterol in persons who combine the use of statins with a high intake of 
phytosterolsZ-stanols, e.g. by the use of functional foods. 

Methods and Results: The proposed model includes the cholesterol pool size in the liver and serum levels of very 
low-density lipoprotein (VLDL) cholesterol. Both an additional and a multiplicative effect of phytosterolZ-stanol 
intake on LDL cholesterol reduction were predicted from the model. The additional effect relates to the decrease 
of dietary cholesterol uptake reduction, the multiplicative effect relates to the decrease in enterohepatic recycling 
efficiency, causing increased cholesterol elimination through bile. From the model, it was demonstrated that a daily 
intake of 2 g phytosterolsZ-stanols reduces LDL cholesterol level by about 8% to 9% on top of the reduction 
resulting from statin use. The additional decrease in LDL cholesterol caused by phytosterolZ-stanol use at the 
recommended level of 2 gZd appeared to be similar or even greater than the decrease achieved by doubling the 
statin dose. 

Conclusion: We proposed a simplified mathematical model to simulate the reduction in LDL cholesterol after 
separate and combined intake of statins and functional foods acting on intestinal (re)absorption of cholesterol or 
bile acids in humans. In future work this model can be extended to include more complex (regulatory) 
mechanisms. 
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Background 

Increased total cholesterol and low-density lipoprotein 
(LDL) cholesterol levels represent a major risk for ather- 
osclerosis and coronary heart disease (CHD). Lipid-low- 
ering drugs, of which the hydroxymethylglutaryl- 
coenzyme A (HMG-CoA) reductase inhibitors (statins) 
have shown to be the most effective, reduce morbidity 
and mortality in patients with CHD [1-3]. Since the last 
decade of the 20"^ century, more interest has been given 



* Correspondence: Cathy.RompelbergSSrivm.nl 

^National Institute for Public Health and the Environment, P.O. Box 1, 3720 
BA Bilthoven, The Netherlands 

Full list of author information is available at the end of the article 



to changing dietary habits, for example with the appear- 
ance of the so-called functional foods. Dairy products 
enriched with phytosterolsZ-stanols are one of the best 
known and most used functional foods to lower elevated 
total and LDL cholesterol levels. PhytosterolsZ-stanols 
are thought to compete with cholesterol for solubilisa- 
tion into mixed micelles, the transport vehicles for cho- 
lesterol across the intestinal wall. Consequently, the 
intestinal (re) absorption of cholesterol is reduced, faecal 
output is increased and total and LDL cholesterol levels 
are lowered by 6% and 10%, respectively [4,5]. Due to 
the rising public awareness of health and nutritional 
improvement, and the mounting evidence of the 
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effectiveness of phytosterolsZ-stanols, it is conceivable 
that in the near future an increasing number of people 
will combine their statin therapy with these functional 
foods. 

To examine the effects on total and LDL cholesterol 
levels of the combined intake of statins and phytoster- 
ols/-stanols, in vivo studies and clinical trials are neces- 
sary. However, for a better interpretation of the 
experimental data as well as to possibly predict choles- 
terol levels given a certain dosing regimen of statins and 
phytosterolsZ-stanols a more theoretically based 
approach is helpful. 

The present study focuses on the combined effect of 
atorvastatin and phytosterolsZ-stanols. However, our 
model can easily be applied to other statins and similar 
acting functional foods (e.g. soluble dietary fibres) as 
well. Moreover, based on certain genetic variants asso- 
ciated with cholesterol absorption and production an 
individual's specific reduction in total and LDL choles- 
terol can be predicted. 

Methods 

We propose a simplified mathematical model to esti- 
mate reductions in LDL cholesterol after separate and 
combined intake of statins and phytosterolsZ-stanols 
(Figure 1). A list of model variables and abbreviations is 
presented in Table 1. Since LDL is the product of very 
low-density lipoprotein (VLDL) delipidation and VLDL 
also transports, although to a lesser extent than endo- 
genous triglycerides, [6] cholesterol from the liver to the 



blood circulation, our model includes the modelling of 
the metabolism of VLDL cholesterol as well. Also a 
hepatic cholesterol pool is accounted for in the model, 
because VLDL cholesterol secretion depends on choles- 
terol pool size. In the next section, we first describe a 
basic model which includes the modelling of the choles- 
terol pool, VLDL cholesterol and LDL cholesterol. Sub- 
sequently, this basic model is reformulated to express 
reductions in LDL cholesterol level dependent of statin 
andZor phytosterolZ-stanol intake. Published scientific 
literature was used to estimate specific model para- 
meters. In the second part of the present study (Results 
section), we tested the appropriateness of our model 
using available published experimental data. 

Basic cholesterol model 
Modelling of the cholesterol pool 

A mass balance is considered with cholesterol input 
from endogenously produced cholesterol Pc and from 
cholesterol taken up from the diet, Ic- Only a fraction 
fabs of dietary cholesterol is assumed to be taken up 
across the gut wall, and consequently the uptake of 
dietary cholesterol is Uc = fabs'^c- The mass balance 
output consists of produced VLDL cholesterol, choles- 
terol cleared by elimination of excess cholesterol 
through bile and cholesterol needed to produce bile 
salts. For simplicity, we neglect the reverse cholesterol 
transport mediated by high-density lipoprotein (HDL) 
and the existence of a hepatic cholesteryl ester pool 
that might be involved in cholesterol homeostasis. 



r \ 

Extrahepatic tissues 



C 



c 



ft 



hack 



V K 

max ' M 



Liver 



Blood 



^ fhack 

-o^ — 



— 

Blood T 




C 

Faeces 

Figure 1 Simplified sclieme of LDL cliolesterol metabolism in humans. For detailed description of the model see text. The definition of the 
model variables are summarised in Table 1. 
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Table 1 Model variables and abbreviations used in the study 


Model variable 


Abbreviation 


EndoQ^nously produced cholcstGrol 


' c 


DiGtBry cho GstGrol intskG 


'C 


ExtcrriBl dsi y stBtin dose 


5 


ExtGrriBl dsi y frGG phytostGrolZ-stsnol dosG 


PS 


{St93dy state) concentration of free cho estero in the liver 


c 


VLDL particles 


V 


Absorbed cholesterol fraction 


fabs 


Fraction of produced VLDL cho estero that re-enters the iver 


'back 


Association rate of VLDL particles and free cho esterol to VLDL cho estero 




Excretion of cholesterol from the cholesterol pool by hi e 




Excretion of cholesterol through the formation of bile salts 


k 

'^salts 


Reduction in cholesterol pool size 


"C 


Cho esterol poo concentration in absence of statins and phytosterolsAstanols 


fn 
^0 


Endogenous produced cholestero in absence of statins and phytosterolsAstanols 


'CO 


Uptake of dietary cho esterol in absence of statins and phytosterolsAstanols 


Urn 


Absorbed cholesterol fraction in absence of statins and phytosterolsAstanols 


f H n 
' abS'O 


Reduction in endogenous cho estero production 


Rn 

np 


Reduction in fraction of cholesterol uptake from the diet 


Rn 


Ratio of exponential rates of different cholestero e imination routes 




VLDL cholesterol production rate 


Pyr 


Transformation rate of VLDL cholesterol to LDL cholesterol 




(Steady state) VLDL cho estero concentration 


VC 


Maximum rate of change in (V)LDL cholesterol due to saturated uptake process 




Michaelis-Menten constant in (V)LDL cholesterol model 


Km 


(Steady state) LDL cholesterol concentration 


LC 


LDL cholesterol production from VLDL cholesterol 


Plc 


Clearance rate of LDL cholesterol through non-saturated process 




Maximal achievable reduction in endogenous cholesterol production 


Pp, max 


Half maximum reduction statin dose 


■5p, 1/2 


Maximal achievable reduction in fraction of cholesterol uptake from the diet 


Pu, max 


Half maximum reduction phytosterolAstanol dose 


P^u, 1/2 



Moreover, up- and down-regulation of LDL receptors 
is not considered. 

The model considers only steady state levels of choles- 
terol, VLDL cholesterol and LDL cholesterol. For a 
steady cholesterol level, the input of cholesterol should 
be balanced by its output. Clearance of cholesterol by 
VLDL cholesterol formation, bile excretion and bile salts 
formation is assumed to be non-saturated and is 
described by the product of clearance rates and steady 
cholesterol level. 

Of the daily amount of VLDL cholesterol formation, 
kvcV -C, the product of steady cholesterol level C with 
VLDL particles V and association rate kyct ^ fraction 
fback is reabsorbed into the liver. It consists of VLDL 
cholesterol that is not used for LDL cholesterol produc- 
tion and of LDL cholesterol. The other fraction 1- fback 
is taken up by the extrahepatic tissues, of which part is 
excreted through HDL cholesterol, which will not be 
considered in this modelling approach. As a 



consequence of this recycling, the effective clearance 
rate of cholesterol to VLDL cholesterol is {'^- fback}'kvc^- 
The amount of cholesterol eliminated through bile salts 
formation is ksaW C. 

Likewise, because it is assumed that not only dietary 
cholesterol but also cholesterol cleared by bile with a 
daily amount of kbUeC is reabsorbed through enterohe- 
patic recycling, the effective clearance rate of cholesterol 
through bile is (I- fabs)-kbiie- 

As we consider the effect of statins and phytosterols/- 
stanols on LDL cholesterol levels, the model becomes 
slightly more complicated. First, it is assumed that 
reduced cholesterol production is related to the 
external daily dose S of statins, Pc = PciS)- Second, it is 
assumed that the reduced cholesterol fraction absorbed 
from dietary cholesterol intake relates to the amount of 
intake of phytosterolsZ-stanols (PS), fabs =fabs{PS)- 

In the steady state the input cholesterol PciS) + fabs 
{PS)-Ic is balanced by cleared cholesterol, which is the 
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product of the effective clearance rates and the steady 
cholesterol \evel({l- fback) kycV) + i^- fabsiPS)) kbiie 
+ksaits)-C{S, PS). Thus, the steady cholesterol level is: 



C{S,PS] ■ 



PciS)+UsiPS)-Ic 



(1 -fback) ■ kvcy + (1 -fabs{PS)) ■ hik + halts 



It should be noted that it is implicitly assumed that 
there is no interaction between statins and phytoster- 
ols/-stanols consumed, i.e. both compounds work simul- 
taneously, independent of each other. 
Modelling of VLDL cholesterol level 

In the modelling of the cholesterol pool (equation (1)) it 
is assumed that the production of VLDL cholesterol Pyc 
is proportional to both the concentration of VLDL parti- 
cles and the free cholesterol level: Pyc = kvcV -CiS, PS). 
Like for cholesterol, a steady state level VC(5', PS) of 
VLDL cholesterol follows from the balance between its 
production and its clearance. VLDL cholesterol is 
assumed to be cleared due to the production of LDL 
cholesterol with daily clearance of kyv^^i^i PS) and 
due to saturated receptor-mediated uptake from blood 
into the liver and extrahepatic tissues [6]. Receptor- 
mediated uptake is assumed to follow Michaelis-Menten 
kinetics with a maximum clearance rate V^ax and a 
saturation constant Km- 

Therefore, a steady state VLDL cholesterol level leads 
to the following mass balance for LDL cholesterol: 



Pvc = kvcV ■ C(S, PS) = kvL ■ VC{S, PS) + 



■ VC{S, PS) 



Km + VC{S. PS) 



Similar as for VLDL cholesterol, a fraction fback of the 
Michaelis-Menten saturated clearance of LDL choles- 
terol from blood goes into the liver. The complementary 
fraction 1- fback goes into extrahepatic tissues (Figure 1). 
(1) ^maxi Km are the same maximum elimination rate and 
saturation constant of the Michaelis-Menten saturated 
uptake of VLDL cholesterol from blood into the liver. 
These constants are assumed to be the same, but this 
assumption is not essential. 

Cholesterol reduction model 
Modelling cholesterol reduction by statins and 
phytosterols/stanols 

When it is assumed that statins and phytosterols/-sta- 
nols work independently of each other, the reduction 
in the cholesterol pool size can be expressed in terms 
of a reduction:7?p (5) = Pc {S)/Pc,o in cholesterol pro- 
duction caused by statin use and a reduction: Ru (PS) 
= fabs{PS)/fabs,o in dietary cholesterol absorption across 
the gut wall caused by phytosterolZ-stanol use. Pc, o 
and fabs,o denote the cholesterol production rate and 
absorption fraction in absence of statins or phytoster- 
ols/-stanols. From equation (1) it is derived in Appen- 
dix 2 that the corresponding reduction Rc{S, PS) in 
the steady state cholesterol level is the product of the 
reduction due to a decrease in enterohepatic efficiency 
(first factor at the right side of equation (4)) and a 
weighted mean of the reduction due to a decrease in 
cholesterol production and uptake (second factor at 
the right side of equation (4)): 



(2) 



The steady state VLDL cholesterol level can be 
obtained by solving the implicit equation for VC(S, PS). 
The explicit expression for VC is deduced in Appendix 
1. Note that of the Michaelis-Menten saturated clear- 
ance of VLDL cholesterol from blood a fraction fback 
goes into the liver. The complementary fraction I- fback 
goes into extrahepatic tissues (Figure 1). 
Modelling of LDL cholesterol level 

In the modelling of the VLDL cholesterol level (equation 
(2)), we assumed that the production of LDL cholesterol 
Pic is proportional to the steady VLDL cholesterol level: 
Plc = kvL-VC{S, PS). LDL cholesterol is assumed to be 
cleared with rate k„ through a non-saturated process 
and by saturated uptake from blood into the liver and 
extrahepatic tissues by the same receptors as for the 
saturated uptake of VLDL cholesterol [6]. Hence, the 
mass balance for steady state LDL cholesterol is: 



Plc = kvLVC(S,PS) = k„ ■ LC{S,PS) + 



V„, 



■LCiS,PS) 



KM+LCiS,PS) 



(3) 



from which the LDL cholesterol level can be obtained 
by solving the implicit equation (3) for LC{S, PS). The 
explicit expression for LC can be found in Appendix 1. 



Rc(S,PS) 



c{s,ps) n + i-fab 



Co 



Ph+1 - fab, (PS) 



Rp(S) 



second tenn \ 

Ru(S) 



1 + Uco/Pco 1+Pco/Uco 



(4) 



second factor 



Here, like for cholesterol production and absorption, 
Co and Uc, o denote, respectively, the cholesterol pool 
concentration and dietary uptake in absence of statins 
and phytosterolsZ-stanols. In the first factor at the right 
side, the ratio denotes the proportion of cholesterol 
elimination through VLDL cholesterol production and 
bile salts production to cholesterol elimination through 
bile excretion, as introduced in Appendix 2. 

The following remarks should be made regarding this 
model. First, the effectiveness of statins or phytosterolsZ- 
stanols to lower cholesterol production is determined by 
the ratio of the contribution of endogenous produced 
cholesterol and the contribution of dietary cholesterol 
uptake to the cholesterol pool. Thus, when dietary cho- 
lesterol uptake is increased, the effectiveness of statins 
(first term in the second factor) is reduced with respect 
to the effectiveness of phytosterolsZ-stanols. Obviously, 
the opposite holds true for the effectiveness of phytos- 
terolsZ-stanols. 
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Second, the reduction in the absorbed fraction of cho- 
lesterol has an additional effect in total cholesterol pool 
reduction (the second term in the second factor at the 
right side of equation (4)) and a multiplicative one (first 
factor at the right side of equation (4)). The additional 
effect relates to the decrease of dietary cholesterol 
uptake reduction, whereas the multiplicative effect 
relates to the decrease in enterohepatic recycling effi- 
ciency, causing increased cholesterol elimination 
through bile. 

Third, the additional reduction caused by statin and 
phystosterolZ-stanol use is a weighted sum of the reduc- 
tion in cholesterol production and cholesterol absorp- 
tion because 1/(1+ Pc, olUc, o) +1/(1+ Uc, o/Pc. o) = 1- 
Formulating reduction as a IVIictiaelis-Menten process 
In the model described above, we aim to associate 
reductions in the cholesterol pool size to reductions in 
LDL cholesterol level. In order to estimate the reduction 
in the cholesterol pool size following statin intake, Rp{S), 
a reduction model has to be assumed. From experimen- 
tal in vitro data from Shum et al [7]. a reasonable 
model assumption is obtained as follows. Shum et al. 
related the concentration of atorvastatin in plasma to 
the inhibition of the enzyme HMG-CoA reductase in 
vitro. We assumed a Michaelis-Menten saturated inhibi- 
tion process. This model was fitted to their experimental 
data and provided a nearly perfect fit (Figure 2). 

From this experimental in vitro result, it is proposed 
that the reduction in endogenously produced cholesterol 
is most likely Michaelis-Menten saturated with adminis- 
tered statin dose as well: 



Rp{S) = 1 



Sp,l/2 + S 



(5) 



where Rp_ ^ux ^ 1 determines the maximum achievable 
reduction and Sp_ 1/2 is the half maximum reduction sta- 
tin dose. 

Based on the fact that cholesterol uptake is receptor- 
mediated, [8] we assumed that the reduction in the frac- 
tion absorbed cholesterol is like: 



J?u(PS) = 1 



■PS 



PSuMi + PS 



(6) 



where R/j, max ^ 1 determines the maximum achiev- 
able reduction and PSu, 1/2 is the half maximum reduc- 
tion phytosterol/-stanol dose. 

Thus, the total reduction in cholesterol pool size after 
combined intake of statins and phytosterols/-stanols is 
obtained by substituting equations (5) and (6) in equa- 
tion (4): 



RclS.PS) -. 



+ 1 - (1 - Ru,^ ■ PS/ (PSu.i/i +PS))- /„te,o 

/, Pc,0 Rp,m«y: S Uc,0 «U,»„ ' K 



Pc,a + Uc,a Sp,i/2 + S Pc,o + Uc,o PSu.i/i + PS 



(7) 



In this derivation we used 1/(1+ Uc, o/Pc, 0)+ 1/(1 + 
Pc. o/Uc. 0) = 1 

Parameter value estimation 

The reduced steady state cholesterol concentration can 
be obtained from a given daily dose of statins or phytos- 
terols/-stanols by applying the reduction model 




Plasma atorvastatin concentration (ng/ml) 

Figure 2 Simulated dose-response relation between in vitro plasma atorvastatin concentration (ng/ml) and hydroxymethylglutaryl- 
coenzyme A (HIVlG-CoA) inhibition (%) in humans. Maximum inhibition is 97.8% and tlie concentration at half maximum inhibition is 1.64 

(ng/ml). Solid symbols present data by Shum et al. [7] 
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proposed in equation (7). From that reduced concentra- 
tion, the reduced steady VLDL cholesterol level can be 
derived by solving equation (2), and subsequently the 
reduced LDL cholesterol level can be derived by solving 
equation (3). However, to be applicable in practice, 
model parameters should be known. In this section all 
parameters of our model are quantified based on data in 
the literature. 

Basic cholesterol model parameters 

It is assumed that the liver produces Pc, o = 1000 mg 
cholesterol per day [9-12]. Furthermore, we assumed a 
dietary cholesterol intake of Ic = 400 mg/ d, of which a 
fraction of 50% (fabs, o = 0.5) [9,13] is taken up in the 
liver (Uc, o = 200 mg/d). The same fraction is supposed 
to be recycled through enterohepatic recycling of cho- 
lesterol excreted with bile. The final contribution to 
liver cholesterol input is assumed to be 70% of produced 
VLDL cholesterol (fback = 0.7) [14]. 

It is assumed that the amount of cholesterol excreted 
with bile is 1000 mg/d[9,13,15] and consequently, 500 
mg/d re-enters the liver. The rate of excretion through 
the formation of bile salts is 400 mg/d [16]. Concerning 
the local liver balance, the input is 1000 (produced cho- 
lesterol, Pc ,o) plus 200 (uptake, Uc, o) plus 500 
(recycled, (I- fabs, o)'kbiie'Co) plus 700 (back transport, 
fback'kvcV'Co) making a total of 2400 mg/d. The output 
is 1000 (bile, kbUeCo) plus 400 (bile salts, ksaitsCo), and 
making a total of 2400 mg/d, plus the production of 
1000 mg VLDL cholesterol per day (VLDL cholesterol, 
kvc^'Co). Moreover, as the elimination from the liver is 
proportional to the production rates of bile salts, choles- 
terol in bile and VLDL cholesterol, the ratio of these 
production rates is ksaits- kbue- kyc^ = 0.4: 1: 1. 

From Sahlin et al. [17,18] we estimated the free cho- 
lesterol content in liver to be 55 nmol/mg microsomal 
protein. Together with a microsomal protein content of 
45 mg/g liver[19] this amounts to 2500 [imol/kg liver 
which equals 960 mg free cholesterol/kg liver. From this 
free cholesterol concentration and the daily bile excre- 
tion, one can derive kbiie = 1000/960 = 1.04, kvcV = 
1.04 and k,aits = 0.416. 

Dietschy et al. [6] report LDL cholesterol model para- 
meter values in humans. When assuming a subject of 70 
kg, these values are V^^x = 1340 mg/d. Km = 90 mg/dl, 
k„ = 5 dl/d and Plc = 910 mg/d. Based on these values, 
a steady cholesterol level LC = 67 mg/dl results from 
equation (3). 

The VLDL:LDL:HDL cholesterol ratio was estimated 
to be 1:8:3 [20]. Thus the corresponding VLDL level is 
8.4 mg/dl. As the LDL cholesterol production rate is 
equal to kyvVC (equation (3)), kvi = 108 (dl/d). From 
equation (2) a VLDL cholesterol production of 1020 
mg/d can be calculated. Above it is assumed to be 1000 
mg/d which shows a consistency error of 2% only. 



Cholesterol reduction model parameters 

From the ratio between the effective clearances, intro- 
duced in the basic cholesterol model parameters section, 
one can derive that the ratio p/^ in equation (4) is 0.7. 
The four remaining parameters Rp_ ^ax, Sp, 1/2, Ru, max, 
PSy, 1/2 are unknown and were fitted to optimise their 
likelihood in comparing modelled LDL cholesterol 
reduction induced by cholesterol reduction to LDL cho- 
lesterol reduction data. Thus, given an estimation of the 
four cholesterol reduction model parameters, the reduc- 
tion in steady state cholesterol is calculated, the result- 
ing reduction in VLDL cholesterol is determined from 
equation (2) and the resulting reduction in LDL choles- 
terol is determined from equation (3). 

To simulate the appropriateness of our model, 
reduced levels are compared with the experimental data 
for separate intakes of atorvastatin[21] and phytoster- 
ols/-stanols [22]. For this procedure we use for Rp^ y„ax> 
Sp^ 1/2 data from a recent meta-analysis of Berry et al. 
[21] that shows experimentally determined in vivo LDL 
cholesterol reduction due to atorvastatin dose. For Ru, 
maxi PSu, 1/2 we use data presented in Demonty et al. 
[22] showing experimentally determined in vivo LDL 
cholesterol reduction due to intake of free phytosterols/- 
stanols, i.e. phytosterols/-stanols not in esterified form. 
Finally, we simulate reductions after combined intake of 
atorvastatin and phytosterols/-stanols. 

Results 

Separate intake 

LDL cholesterol reduction by atorvastatin 

We applied equation (7) together with the correspond- 
ing VLDL and LDL cholesterol levels equations (2) and 
(3) to data in Berry et al. [21] in a LDL cholesterol 
reduction model using the model parameters given 
above. The unknown parameter values in equation (7) 
were estimated through fitting the maximum reduction 
Rp, max aud the statin dose Sp^ 1/2 when half maximum 
reduction is reached. 

Figure 3 shows the comparison of the resulting LDL 
cholesterol reduction model to the data in Berry et al. 
[21] Parameter values are Rp^ ^ax = 0.544 (standard error 
(SE) = 0.033) and half reduction dose Sp^ 1/2 = 6.7 (SB = 
1.4) mg/d. Given the parsimony of dose levels and scat- 
tering of data, a good comparison between the model 
and the experimental results is obtained (R^ = 0.70). 
LDL cholesterol reduction by phytosterols/stanols 
We applied equation (7) together with the correspond- 
ing VLDL and LDL cholesterol levels equations (2) and 
(3) to data in Demonty et al. [22] in a LDL cholesterol 
reduction model. The unknown parameter values in 
equation (7) were estimated through fitting the maxi- 
mum reduction 7?^ ^^a: and the free phytosterol/-stanol 
dose PSu 1/2 when half maximum reduction is reached. 
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Atorvastatin dose (mg/d) 

Figure 3 Simulated reduction (%) in LDL cholesterol after treatment with different doses of atorvastatin The solid line sliows tine fit to 
the model and symbols represent experimental data from Berry et al. [21] Values for the Michaelis-Menten parameters are: effective maximum 
LDL cholesterol reduction (ftp^ ^^x) = 0.544 and half maximum reduction statin dose (S^ ,/2) = 6.7 mg/d. The dashed lines show the 5% and 95% 
uncertainty range in reduction obtained by correlated sampling (correlation coefficient p = 0.88) of fip, ,„ax and Sp v2 froni their covariance 
matrix. 



Figure 4 shows the comparison of the resulting LDL 
cholesterol reduction model to the data of Demonty et 
al. [22] Parameter values are Ru, max = 0.221 (SE = 
0.039) and half reduction dose PSu, 1/2 = 1-78 (SE = 
0.69) mg/d. The model shows reasonable agreement 
with the published experimental data (R^ = 0.17). 



Combined intake 

LDL cholesterol reduction by combined use of atorvastatin 
and phytosterols/stanols 

The model is applied to LDL cholesterol reduction due 
to the combined intake of statins and phytosterols/-sta- 
nols. For subjects with a daily intake of 0, 20, 40 and 80 




PhytosterolZ-stanol dose (g/d) 

Figure 4 Simulated reduction (%) in LDL cholesterol after treatment with different doses of phytosterolsZ-stanols. The solid line shows 
the fit to the model and symbols represent experimental data from Demonty et al. [22] Values for the Michaelis-Menten parameters are: 
effective maximum LDL cholesterol reduction (Ru, max) = 0.221 and half maximum reduction phytosterolZ-stanol dose {PSu, 1/2) = 1-78 mg/d. The 
dashed lines show the 5% and 95% uncertainty range in reduction obtained by correlated sampling (correlation coefficient p = 0.98) of Ru, max 
and PSu, 1/2 from their covariance matrix. 
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mg atorvastatin, respectively, we show in Figure 5 the 
total reduction in LDL cholesterol as a function of daily 
phytosterolZ-stanol intake. For the daily recommended 
intake level of 2 g free phytosterolsZ-stanols (equivalent 
to 3.3 g/d phytosterol/'Stanol esters), the additional 
reduction is 4.2% for a daily statin dose of 80 mg, 4.5% 
for a daily statin dose of 40 mg, 4.8% for a daily statin 
dose of 20 mg and 7.8% with no statin intake. Thus, the 
reduction in LDL cholesterol caused by additional phy- 
tosterolZ-stanol intake decreases with increasing daily 
atorvastatin dose. However, when considering the addi- 
tional decrease as a percentage of the LDL cholesterol 
level already reduced due to statin intake, the additional 
decrease ranges from 7.8% (no statin intake) to 8.6% (80 
mg daily statin dose). Randomised controlled trials in 
which patients on statin therapy were treated daily with 
1.8 to 6 g phytosterol or -stanol esters have shown 
reductions in the same order of magnitude, i.e. between 
6.1% and 10.3% [23-28]. 

The reduction in enterohepatic recycling contributes 
for 68%, 58%, 56% and 55% of the total decrease in 
LDL cholesterol levels at daily statin doses of 0, 20, 40 
and 80 mg, respectively. At the same recommended 
phytosterolZ-stanol intake level of 2 g/d, the additional 
decrease in LDL cholesterol by physterolsZ-stanols for a 
daily statin dose of 20 mg (4.8%) is equal to the addi- 
tional decrease by doubling daily statin dose to 40 mg 
(5.3%). For a daily statin dose of 40 mg the additional 
decrease in LDL cholesterol by physterolsZ-stanols 
(4.5%) is 30% larger than the additional decrease by 
doubling daily statin dose to 80 mg (3.2%). 



Discussion 

In this paper, a mathematical model is presented that 
simulates the reductions in LDL cholesterol after sepa- 
rate and combined intake of atorvastatin and phytoster- 
ols/-stanols in humans. We demonstrated that a daily 
intake of 2 g phytosterolsZ-stanols reduces LDL choles- 
terol level by about 8% to 9% on top of the reduction 
resulting from statin use. This level of reduction is con- 
sistent with the findings of randomised controlled trials 
[23-28]. The additional decrease in LDL cholesterol 
caused by phytosterolZ-stanol use at the recommended 
level of 2 gZd appeared to be similar or even greater 
than the decrease achieved by doubling the statin dose, 
a finding that has been observed previously in human 
trials [23,29]. The reduction in LDL cholesterol level 
due to phytosterolZ-stanol use results from a decrease in 
the intestinal uptake of dietary cholesterol (additional 
effect) and a reduction in enterohepatic recycling (multi- 
plicative effect). For daily statin doses of 20 mg or more, 
the contribution of the enterohepatic recycling reduc- 
tion is 55% or more. When no statin is used, this contri- 
bution is 68%. 

Mathematical models provide a valuable means of 
interpreting experimental data and improving the ability 
to predict the response to a given treatment. Other 
modelling studies have focused on cholesterol metabo- 
lism, but are merely aimed at answering questions on 
the cellular level or tend to focus on specific areas of 
cholesterol metabolism, such as the fluid dynamics of 
lipid accumulation on the arterial wall or the chemical 
kinetics of LDL oxidation [30-32]. 
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In the present study, the separate and combined 
effects of the cholesterol-lowering drug atorvastatin and 
functional foods with phytosterolsZ-stanols in humans 
were modelled. Yet, this model can easily be applied to 
other statins and similar acting (functional) foods as 
well. Products with soluble dietary fibres, for example, 
are also known to lower total and LDL cholesterol by 
reducing the intestinal (re)absorption of cholesterol and 
bile acids, although they work by a different mechanism 
as phytosterolsZ-stanols [4,33,34]. Moreover, individuals' 
specific reductions in total and LDL cholesterol can be 
predicted, based on certain genetic variants. For exam- 
ple, the ratio of cholesterol synthesis to cholesterol 
absorption varies between individuals and is an impor- 
tant determinant for the cholesterol pool size [35]. Also 
mutations in the LDL receptor gene causing familial 
hypercholesterolaemia can be modelled by varying the 
parameter V^„. 

There are a few possible directions for improving our 
model. First, the model could be extended by including 
the up- and down regulatory mechanisms involving the 
LDL receptors. Nonetheless, since we assumed that the 
clearing of {V)LDL cholesterol from the blood follows 
Michaelis-Menten kinetics, we implicitly included recep- 
tor-mediated uptake in the model. Also other regulatory 
control pathways were disregarded, such as the existence 
of a hepatic cholesteryl ester pool that might be 
involved in cholesterol homeostasis and the regulatory 
loop in the synthesis of LDL receptors [30]. Another 
extension would include reverse cholesterol transport 
mediated by HDL [20]. Moreover, the proposed model 
assumes that the reducing effects of statins and func- 
tional foods are independent of each other. Although 
this is likely the case for the combination of phytoster- 
olsZ-stanols and statins, [22,24,36,37] it is uncertain 
whether this applies for other food-drug combinations. 
It has, for example, been proposed that soluble dietary 
fibres reduce the intestinal uptake of statins [38,39]. Our 
model should be extended to include such an 
interaction. 

In conclusion, we proposed a simplified mathematical 
model to simulate the reduction in LDL cholesterol 
after separate and combined intake of statins and func- 
tional foods acting on intestinal (re)absorption of cho- 
lesterol or bile acids in humans. In future work, this 
model can be extended to include more complex (regu- 
latory) mechanisms. 

Appendix 1 

Steady VLDL cholesterol concentration 

The mass balance equation (2) can be rewritten as a 
quadratic equation in VC: 



fevL ■ VC^ + (V„ + IzvL ■ Km - kvcV QVC-Km- kvcV 0 = 0 

Here, we suppress the dependencies on statin admin- 
istration S and dietary phytosterolZ-stanol intake PS. 

A quadratic equation has two solutions, but the only 
physicochemical relevant solution for which the VLDL 
cholesterol concentration is non-negative is: 

VC = -|- {kvcV ■ C - kvL ■ Ku - V™ 

+J (kvcV -C-kvL-Ku- V„^f + 4kvL ■ Km ■ kycV ■ c] 



Steady LDL cholesterol concentration 

Like for VLDL cholesterol, the mass balance equation 
(3) can be rewritten as: 

k„ ■ LC^ + (V,„,„ + fe„ • Km - kvL ■ VC] ■LC-KM-kvL-VC = 0 
with as solution: 

I-C = (fevL -VC-K-Km- Vmax 



+V (kvL -VC-kn-Ku- Vmaxf + 4fe„ ■ Km ■ kvL ■ VC 



Appendix 2 

Reduction in steady state cholesterol 

The reduction in steady state cholesterol level is: 



^ {Pels) + /;,h(PS) - Jc) / (( 1 - fb^t] ■ kvcV + (1 - /a(PS]] ■ hii, + 

{Pc,0 +/abs,0 ■ Jc)/{(1 -fiack) ■ ^VcV + (1 - fahs,o) ■ hih + ^salt) 

^ (1 -had,) ■ ^tycV * (1 -/a,o) ■ hih * fe«il8 Pc(S) +/a(PS) ■ Jc 

(1 -had,) ■ fevcl'+ (1 -/a(PS)) ■ few, + kMio Pc,o +/afa,0 ■ 'c 
^ Pt+l-/A.o / Pc(S) ^ U,[PS) lc \ 

Pk + l-hbs{PS) \Pc.O + Uc,0 PC.O + hbs.O ■ tc / 

Pt+l-/A.o / Rp(S) ^ Ru(PS) \ 
» + 1 - /a(PS) a 1 + L7c,o/Pc,o * 1 + Pqo/Uqo ) 

In the third line, the ratio of clearance rates = ((1- 
fback)-kvcV + ksaits)lkbiie is introduced. In this line also 
one instance oifabs, o'lc is substituted by Uc, o- In the 
fourth line the definition of production reduction and, 
after dividing out intake Ic, of uptake reduction is 
substituted. 
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